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Abstract 

A recent study of coherent noise model for the system size independent case pro- 
vides an exact relation between the exponent r of avalanche size distribution and 
the q value of the appropriate g-Gaussian that fits the return distribution of the 
model. This relation is applied to Ehrenfest's historical dog-flea model by treating 
the fluctuations around the thermal equilibrium as avalanches. We provide a clear 
numerical evidence that the relation between the exponent r of fluctuation length 
distribution and the q value of the appropriate g-Gaussian obeys this exact rela- 
tion when the system size is large enough. This allows us to determine the value 
of g-parameter a priori from one of the well known exponents of such dynamical 
systems. Furthermore, it is shown that the return distribution in dog-flea model 
gradually approaches to g-Gaussian as the system size increases and this tendency 
can be analyzed by a well defined analytical expression. 
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1 Introduction 



Since tlie pioneering work of Bak, Tang, and Wiesenfeld (BTW) in 1987 [T], 
the self-organized criticality (SOC) paradigm has seen a burst of activity in the 
literature. In their seminal paper, by making use of so-called BTW sandpile 
model, these authors demonstrated that systems under their natural evolution 
are driven at a very slow rate until one of their elements reaches a critical 
stationary state and this triggers a sudden activity, i.e., avalanche. When the 
avalanche is over, the system evolves again according to the slow drive until a 
next avalanche is triggered without having any fine-tuning of parameters. Such 
a dynamic gives rise to the power-law correlations seen in the non-equilibrium 
steady states [2|3] . 

Following the BTW sandpile model a great variety of models ranging from 
the deterministic and stochastic to the dissipative and conservative ones have 
been introduced which exhibit the phenomenon of SOC (for an overview, see 
[1] and references therein). Among them, a random neighbor version of the 
BTW sandpile model shows that a dynamical system with only two degrees 
of freedom can exhibit SOC and the dynamics can be described by a master 
equation [5]. Soon after the random neighbor BTW sandpile model was in- 
troduced, its conservative variant has been reported piT] . This conservative 
variant is neither extended nor dissipative with regard to the amount of sand 
in the system but still shows SOC with critical exponents where the dynamics 
of the model is given by a Fokker-Planck equation. The avalanche size distri- 
bution P(A) is readily obtained by solving the Fokker-Planck equation at an 
absorbing boundary and is shown to exhibit a power-law regime -P(A) ~ A""^, 
followed by an exponential tail. Indeed, this model is an adaptation of the 
famous dog-flea model introduced by Ehrenfest in 1907 [S] to describe the 
process of approaching an equilibrium state in a large set of uncoupled two 
state systems in the presence of fluctuations (avalanches) around this state 

In a recent work, the SOC feature of the dog-flea model was studied by sim- 
ulating the underlying stochastic process that describes the natural evolution 
of the model [Tlj. In this paper, we show that the relation between the power- 
law exponent r of the fluctuation length distribution and the q value of the 
appropriate g-Gaussian that fits the return distribution {i.e., distribution of 
fluctuation length differences at subsequent time steps) obeys the rule 



which was reported for the limited number of earthquakes from the World 
and California catalogs [12]. This approximate relation ([1]) enables one to 
obtain the value of q parameter a priori from the power-law exponent r of 
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fluctuation length distribution. Using the same line of thought, the r value for 
the dog-flea model was obtained as r ^ 1.517 by making use of the maximum 
likelihood estimation method yielding to a value of g = 2.35 of the concomitant 
g-Gaussian that fits the return distribution at iV — )• oo limit [TT] . 

Soon after the relation ([1]) was reported, an exact relation between r and q has 
been introduced for the coherent noise model (CNM) in the size independent 

case as 



(2) 



These relations between r and q given by Eqs. ([T]) and ([2]) are slightly different 
from one another. Both relations approach g = 1 at r — ?■ oo limit while only 
the latter achieves the value g = 3 for the given r = 1 value [13] . Let us 
also remark that if we identify r = l/(g5 — 1), we can obtain from Eq. ([2]) 
gs — 1 = (g — l)/2, which is precisely Eq. (30) in Ref. [I] 



The chosen system size plays a crucial role in analyzing the SOC feature of 
the model. At relatively small system sizes, the distribution of the fluctuation 
length time-series \{t) in dog- flea model exhibits a power-law regime, which 
is followed by an exponential decay because of the finite size effects, while at 
large system sizes, the power-law regime increases and the exponential decay 
is postponed [11]. This size dependent behavior of the fluctuation length dis- 
tributions results in different topological properties of the return distributions 
as we discuss in subsequent sections. 

Our main task will be to analyze the SOC in the dog-flea model through 
numerical evaluation of the underlying stochastic process. In this manner the 
fluctuation length and return distributions are studied when the system size 
is large enough, i.e., at the N oc limit, and for small system sizes {i.e., 
when the finite size effects are visible). The rest of the paper is organized as 
follows. The model and the numerical procedure that we implement are given 
in the next section. Then, the probability distribution of fluctuation length for 
the system size, chosen so that the finite size effects are avoided, is obtained 
by numerical evaluations. Once the power-law exponent r of the fluctuation 
length distribution is found, the distribution of returns are analyzed. Next, 
we analyze the return distribution of the dog-flea model when the finite size 
effects are visible. A summary and discussion of the results conclude this 
communication. 
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2 Dog-flea model 



The dynamics of the dog-flea model has simple rules. The model has dy- 
namical sites represented by the total number of fleas shared by two dogs (dog 
A and dog B). Suppose that there are Na fleas on dog A and Nb fleas on dog B 
leading to a population of fleas N = + Nb- For convenience, N is assumed 
to be even. In every time step, a randomly chosen flea jumps from one dog to 
the other. Thus, we have Na — )• A^^ ± 1 and Nb — >■ Nb =F 1- The procedure 
is repeated for an arbitrary number of times. In the long time run, the mean 
number of fleas on both dog A and dog B converges to the equilibrium value, 
(Na) = (Nb) = N/2 with the fluctuations around these mean values. A single 
fluctuation is described as a process that starts once the number of fleas on 
one of the dogs becomes larger (or smaller) than the equilibrium value N/2 
and stops when it gets back to the value N/2 for the first time. Thus, the end 
of one fluctuation specifies the start of the subsequent one. The length (A) 
of a fluctuation is determined by the number of time steps elapsed until the 
fluctuation ends. 



3 Fluctuation length and return distributions 

In order to study the behavior of return distributions with large A^ the system 
size is chosen as A^ = 10®, enabling the finite size effects to appear only 
at very large fluctuation lengths. In Fig. [TJ^a), we plot the distribution of 
the fluctuation length time-series A(t) obtained from 2 x 10^ fluctuations. It 
can then be seen that the fluctuation length distribution follows a power-law 
regime, P{X) ~ A"^ with r ~ 1.517 and exponential decay {i.e., finite size 
effects) is not visible. 

In Fig. dl^b), the distribution of returns, i.e., the difference between fluctuation 
lengths obtained at consecutive time steps, as AA(t) = A(t + 1) — A(t) is shown. 
At this point it should be emphasized that in order to have zero mean, the 
returns are normalized by introducing the variable a; as a; = AA — (AA), where 
(■ ■ ■) stands for the mean value of the given data set. From Fig.[T](b), it is clear 
that the signal of return distribution for A^ = 10® is not Gaussian. Instead, it 
exhibits fat tails which is described by g-Gaussian, 

P{x) = P(0)[(1 + Piq - l)a;2]V^''^\ (3) 

where /3 characterizes the width of the distribution and q is the index of 
nonextensive statistical mechanics [T51IT6] . As suggested in [13], the parameter 
q can be determined directly from the relation ([2]) a priori as q = 2.32. In 



4 



lo' 10" 10^ lO'* 

A, 




-500 -250 250 



xP(0) 

Fig. 1. (Color online) (a) Fluctuation length distribution for = 10^. The full 
black line represents the fitting curve of the distribution with slope r ~ 1.517. The 
distribution has an arbitrary normalization such that P{X = 1) = 1. (b) Return 
distribution for N = 10^. The full black curve represents the appropriate g-Gaussian 
with q = 2.32 determined from relation ([2]) a priori and /? = 30. The central part 
of the distribution is emphasized in the inset. 

Fig. mb), the solid black line represents the appropriate g-Gaussian with q = 
2.32 and /3 = 30 that perfectly fits the signal of return distribution not only 
for the tails but for the intermediate and the very central part (see inset). 

The fluctuation length distribution P{X) in dog-flea model shows size depen- 
dent behavior. Its signal exhibits a power-law regime following by an expo- 
nential regime because of finite-size effect at relatively small system sizes. 
However, as it has been reported in [TT], the finite-size scaling hypothesis is 
satisfied. In Fig. [2]^a), the fluctuation length distribution of 10^ fluctuations 
for four different system sizes N = 10^, N = 10^, = 5 x 10^ and N = 10"^ 
are presented. As the system size increases the power-law regime gets longer 
and the exponential decay is postponed. Since the system size is far away 
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Fig. 2. (Color online) (a) Fluctuation length distributions for = 10^, = 10'^, 
AT = 5 X 10^ and iV = 10^. 10^ fluctuations are considered for each different system 
sizes. The legends of figure give the corresponding slope values r of power-law 
regimes. The distributions have an arbitrary normalization such that P{X = 1) = 1. 
(b) Return distributions for iV = 10^, iV = 10^, iV = 5 x 10^ and = 10^. The black 
lines on each return distribution represent the fitting curves obtained by Eq. (j4j) with 
q values determined a priori from the relation ^ using the corresponding r values. 
For each curve ag = 30 which is equal to /3 value of appropriate g- Gaussian that 
fits the return distribution with largest system size, i.e., N =10^. In the inset, we 
also give the parameter oi as a function of A^ in order to better visualize how the 
finite-size effect disappears as A^ increases. It is clear that ai appears to vanish as 
ai = 963A^-2_ 

from being at the N ^ oo limit {i.e., not large enough to avoid the finite 
size effects), each distribution of fluctuation length has different power-law 
slopes with slightly different r values. This difference between the r values 
of each distribution of fluctuations for different system sizes leads to return 
distributions, which are neither Gaussian nor g-Gaussian. Nevertheless, from 
Fig. Mjo) it is obvious that as the system size N increases, leading to a longer 
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power-law regime in the fluctuation length distribution, the curves of return 
distributions start to exhibit a tendency to a kind of fat tailed distribution, i.e., 
g-Gaussian. This gradual approach to g-Gaussian in the presence of finite-size 
effect is described by the following distribution 



a. 



ai ai 



'qjq-l)aix'^ 



1/(1-9) 



(4) 



If ai = 0, Eq. (jl]) coincides with the g-Gaussian and a Gaussian is recovered 
for Gg = ai {i.e., q = 1). 

As it is the case in CNM [13], Eq. (jl]) seems to coincide with the signal of 
the return distribution when the finite-size effects are invisible in dog-flea 
model. As it is seen from Fig. [2t^b), each return distribution corresponding 
to different system sizes exhibits a topological behavior which is fitted by 
Eq. 01]) with appropriate q values determined a priori from the relation ([2]). 
Each critical exponent r for different system sizes is obtained by considering 
only the power-law regime of corresponding fluctuation length distribution, 
see Fig. [2]^a). Fig. Mjo) also reveals the fact that the longer the power-law 
regime persists for fluctuation length distribution, the better the appropriate 
g-Gaussian dominates in the return distribution. Eventually, as — )■ oo, the 
power-law regime of fluctuation length distribution is expected to continue 
forever (see Fig. [TJ^a)). It is also evident from the inset of Fig. Et^b) that the 
finite-size effects vanish as A^~^, which is completely the same tendency found 
in [13j for the coherent noise model. When the finite-size effects disappear, the 
corresponding return distribution seems to converge to the g-Gaussian for the 
entire region (see Fig. [Dl^b)). 



4 Conclusion 



We analyze the SOC in Ehrenfest's dog-flea model through the probability dis- 
tributions of fluctuation length and of the differences between the fluctuation 
lengths at subsequent time steps, i.e., returns, by simulating the stochastic 
process that describes the evolution of model. The fluctuations around the 
thermal equilibrium are treated as avalanches. In order to avoid the finite- 
size effects the size of the system is chosen large enough, i.e., N = 10^. This 
enables one to obtain a power-law regime with slope r ^ 1.517 without any 
exponential decay. Then, the signal of return distribution is analyzed and it 
is shown that it converges to a g-Gaussian with g = 2.32, a value obtained a 
priori from Eq. (|2]). This g value is slightly different from the one reported in 
pT| . where g = 2.35 has been obtained using Eq. ([1]). The difference between 
the g values for the same distribution stems from the fact that Eq. ([T]) is an 
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approximate relation slightly differing from the exact relation given by Eq. ([2]) 
(see [13] for details). 

The case where the finite-size effects are visible is also investigated by exten- 
sive simulations. When the system size is chosen relatively small, the finite- 
size effects appear in the avalanche size distribution: the distribution follows a 
power-law regime with corresponding slope value r followed by an exponential 
decay. In this case, the return distributions numerically converge to appropri- 
ate g-Gaussian starting from the central part and gradually evolves towards 
the tails as the system size increases. This behavior of return distribution is 
indeed in good agreement with the one obtained when the system size is large 
enough. It is also found that the finite-size effects vanish as A^~^. As — )■ oo 
the signal of return distributions for the finite system sizes numerically con- 
verges to a g-Gaussian in the entire region. 
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